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Abstract. We present a complete survey of current star formation in the Perseus molecular cloud, made at 850 and 450 pm 
with SCUBA at the JC MT. Covering 3 deg 2 , th is submillimetre continuum survey for protostellar activity is second in size only 
to that of p Ophiuchus l Johnston e et all2004l) . Complete above 0.4 M Q (5<x detection in a 14" beam), we detect a total of 91 
protostars and prestellar cores. Of these, 80% lie in clusters, representative of star formation across the Galaxy. Two of the 
| groups of cores are associated with the young stellar clusters IC348 and NGC1333, and are consistent with a steady or reduced 

■^j- , star formation rate in the last 0.5 Myr, but not an increasing one. In Perseus, 40-60% of cores are in small clusters (< 50 M G ) 

l/^ ' an d isolated objects, much more than the 10% suggested from infrared studies. Complementing the dust continuum, we present 

a C ls O map of the whole cloud at 1' resolution. The gas and dust show filamentary structure of the dense gas on large and small 
scales, with the high column density filaments breaking up into clusters of cores. The filament mass per unit length is 5-1 1 M Q 
per 0. 1 pc. Given these filament masses, there is no requirement for substantial large scale flows along or onto the filaments in 
order to gather sufficient material for star formation. We find that the probability of finding a submillimetre core is a strongly 
O f increasing function of column density, as measured by C 18 integrated intensity, P(core) oc I 3 , This power law relation holds 

down to low column density, suggesting that there is no A v threshold for star formation in Perseus, unless all the low-A v submm 
cores can be demonstrated to be older protostars which have begun to lose their natal molecular cloud. 
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1 . Introduction the point of formation or shortly after, at protostellar cores and 

5h ' deeply embedded Class objects, with ages less than or of or- 

Cd Recent advances in submillimetre (submm) detector technol- der 1q4 years wkh submm suryeys we can rdate the positions 

ogy now make it possible to image entire star forming regions of ±e young protostars to the p r0 p er ties of the dust and gas 

and detect all the protostellar and starless cores within. For the in the m olecular clou ds in which they form. This has led other 

first time, we can gather statistically significant samples of stars amhors J 0nish j eta lJ|l99Rt I.Tohnstone et aljEool to conclude 

at the earliest stages of their evolution. At the same time, theo- ^ there . g an extinction or column density thresh old for star 

nes of star formation are moving on from the formation of indi- formation? that ^ that stars only form m me densest parts of 

vidual, isolated ejects (e.g.,|Shu et al.||1987) to the formation molecular clouds We investigate this conc i usion furth er in this 

of clusters (e.g.. lMac Low & Klessenll2004l) . It is thus hugely p a p er 

important at this time to make large surveys of star formation 

regions to compare with the theories. We have carried out a complete survey of star formation 

An important question which we can only address with above 0.4 M Q in the Perseus molecular cloud. Perseus was 

submm surveys is where in molecular clouds stars form. This selected as it is known to be forming clusters of stars and 

cannot be answered by looking at older populations of pre- more massive stars than, e.g., Taurus. The Perseus molecular 



main-sequence stars detectable in the infrared and optical be- cloud contains a total of 17,000 M n of gas, es timated from 
cause the typical velocities of the stars lead them to wander 



cloud contains a total ot 1/,UUU M n ot gas, estimated from 
visual extinction jBachiller & Cernicharolll986bl) . Our survey 
by 1 pc/Myr, and the environment in which they formed may area of ~ 3 deg 2 was selected to include all the moderate 
no longer exist. Therefore we are driven to look for stars at to high extinction material (A v > 4) where star formation 

might occur. Imaging dust emission at 850 and 450 pm with 

Send offprint requests to: hatchell@astro.ex.ac.uk the Submillimetre Common User Bolometer Array (SCUBA) 
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at the James Clerk Maxwell Telescope (JCMT), we detect 
all the protostars and starless cores above 0.4 M Q /14" beam. 
Complementing the dust continuum images, we present molec- 
ular line (C 18 7=1-0) data at 1' resolution. Thus, we can re- 
late the star formation activity to the underlying physical prop- 
erties of the molecular cloud and to the kinematics. By imaging 
a single molecular cloud we minimise the uncertainties in de- 
rived masses, sizes, etc., due to distance. 

The Perseus molecular cloud is a well known star forming 
cloud in the nearby Galaxy. It is associated with two clusters 
containing pre-main-sequence stars: IC3 48, with an estimate d 
age of 2 Myr (and spread of ±1.5 Mvr: Luhman et alJl2003l): 
NGC1333, which is less than 1 Myr in age JLada et al.lll996l 
IWilking et al]l2004|): and the Per 0B2 asso ciation, which con- 
tains a BO. 5 star Jsteenbrugge et al^OO ^ and therefo re must 
be less than 13 Myr in age jMevnet & Maederll200ol) . These 
clusters therefore all show evidence of star formation activity 
within the last ~ 10 7 years. The molecular cloud itself contains 
a number of previously known protoclusters and isolated pro- 
tostars but until now there has been no complete census of the 
current star formation activity. 

We assume the distance to the Perseus molecular cloud to 
be the same as the Hipparco s distance to the embed ded clusters, 
and use 320 pc throughout dde Zeeuw etalJl999t) . 



2. Observations and data reduction 

Images of Perseus were obtained at 850 and 450 //m us- 
ing SCUBA on the JCMT during 20 nights between 1999 
and 2003. The region selected corresponds roughly to the 
fT* A dV = 4Kkms~ 1 contour of the Bell labs 7m 13 CO map 1 , 
although not all of this area was mapped with both SCUBA and 
C ls O, while some regions outside this contour were mapped 
with each. For completeness and calibration consistency we 
included the NGC1333 region although i t was m apped pre- 
viously with SCUBA bv lSandell & Kneel ( EoOlal) . The whole 
region covered contains more than 10,000 M of gas, esti- 
mated from 13 CO. Fields of size 10' x 10' were scan-mapped 
6 times with each of 3 chop throws (30", 44" and 68") and 
2 chop dir ections (R.A. and Dec.). Th e data were reduced us- 
ing SURF Jjenness & Lightfoolll99'8l The typical RMS noise 
level is 35 mJy/beam at 850 /urn. The weather conditions for 
the 450 yum observations were more variable: for the clus- 
ters (observed in the best conditions as we expected to find 
sources here) the noise level is ~ 200 mJy/beam, but increases 
to > 1 Jy/beam in other regions. The 450 and 850 /jm beam 
sizes are 8" and 14" respectively or 0.012 and 0.022 pc at the 
distance of Perseus. We assum e a dust emissivity of / ^so = 
0.012 cm 2 g -1 (gas+dust), from lOssenkopf & Henningl l l 19941) 
model 5 for icy coagulated grains and a gas-to -dust mass ra- 
tio of 161 (1.5 x 10~ 26 g of dust per H at om: iDraine & Lee 
1984 ), assuming gas is 0.89% H by number JChurchwell et al. 
19901) . Our beam mass sensitivity to a 12 K core is then 0.4 M Q 
at 850 yum, although we also detect warmer, less massive ob- 



jects, and our beam-averaged column density sensitivity is 
6 x 10 22 cirT 2 (both limits given for 5cr detections). 

Maps of the molecular cloud in C ls O /=l-0 were made in 
January and December 2000 using the 4x4 array SEQUOIA 
on the FCRAO 14m telescope. The beamsize was 46", corre- 
sponding to 0.07 pc at the distance of Perseus. The 92,500 spec- 
tra were baselined, calibrated to 7\ib with an efficiency of 0.45, 
smoothed to 0.25 km s _1 resolution, gridded and integrated in 
velocity over 0-12 km s _1 using the GILDAS software, to pro- 
duce the integrated intensity map shown in Fig.[0 Visualisation 
and further analysis of th e C 18 Q were carried out using the 
Karma visualisation suite JGoochlll996l ). Maps of 13 CO(1-0) 
were also obtained during this same observing period, the re- 
sults from which will be presented at a later date. 

We take jV(C 18 Q) = 3.4 x 10 I4 (A V - 1.6) cirT 2 

dBachiller & Cernicharolll986bl) . or alternatively, 7(C 18 0) = 
(A v - 1.6)/3.6 K km s" 1 , assuming LTE at 12 K. We also 
assume 7Y(H+2H 2 ) = 0.9 x 10 20 A V dBohlin etalJll978h . The 
temperature of the molecular gas is typically 10— 12 K from 
NH3 observations dBachiller & Cernicharolll986al). although 
protostell ar cores are warmer (30-50 K: lMotte & Andrel2"oOU 
Ijennings et al.ll987l) . 

Our definition of a core for the purposes of this paper is 
a flux peak with closed contours at the 500, 250 or 150 mJy 
(~ 4<x) levels. Elongated regions at those levels were divided 
into subpeaks at subsidiary maxima. It should be borne in mind 
that it is not necessarily straightforward to convert our obser- 
vational sensitivity of 0.4 M /beam to the minimum total mass 
of a core identifiable in our survey as this limit depends on the 
core density distribution and size: massive but extended cores 
could easily remain undetected. As a consequence it is likely 
we are least sensitive to the youngest and hence least centrally 
condensed prestellar cores. In addition if the pressure external 
to a core determines its size for a given mass, our core mass 
sensitivity may vary from region to region in the cloud giv- 
ing a hi gher limit in the lower pressure regions away from the 
clusters (iTohnstone et al.l2004l) . We are also insensitive to weak 
sources in the emission wings of strong ones and therefore un- 
derestimate the number of weak sources in clusters, an effect 
which is discussed further in the next section. 



1 The 13 CO map is available from John Bally, at 
http : //casa . Colorado . EDU/~bally/f tp/7mdata/ 
Perseus_13co. fits 
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Fig. 1. a. The eastern part of the Perseus molecular cloud complex. Top: SCUBA 850 //m map. The colourscale is from to 
2281 mJy/beam on a square root (flux density) scale to bring out the low intensity features. Bottom: Contours of C 18 1-0 
integrated intensity at 1' resolution overlaid on 850 fim colourscale. Contours are 1, 2, 4, 8, 16 K km s _1 . In both maps the black 
boundary marks the area mapped with SCUBA and the grey boundary marks the area mapped in C ls O (note, B5 was not mapped 
in C 18 0). Some well-known clusters and IRAS sources are marked. 
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Fig. 1. b. The western part of the Perseus molecular cloud complex. Top: SCUBA 850 fim map. The colourscale is from to 
2281 mJy/beam on a square root (flux density) scale to bring out the low intensity features. Bottom: Contours of C 18 1-0 
integrated intensity at 1' resolution overlaid on 850 /im colourscale. Contours are 1,2, 4, 8, 16 K km s" 1 . In both maps the black 
boundary marks the area mapped with SCUBA and the grey boundary marks the area mapped in C ls O (note, B5 was not mapped 
in C 18 0). Some well-known clusters and IRAS sources are marked. 
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Fig. 2. SCUBA 850 jum maps of the (a) B5, (b) HH21 1, (c) Bl, (d) L1448, (e) NGC1333, (f) L1455, (g) IRAS 03292+3039, (h) IRAS 03282+3035 and (i) IRAS 03271+3213 
regions. Colourscale is 850 fim continuum from (white) to 2281 mJy/beam (black) on a square root (flux density) scale to bring out the low intensity features. Black contours 
are 850 /im 100, 200, 400, 800, 1600, 3200 mJy/beam. White contours are 450 yum 2, 4, 8 Jy/beam. Dark green contours are where the C 18 integrated intensity = 1 K km s _1 . 
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3. Clusters 

The complete 850 /an image in Fig. n snows emission from 
dust cores and high column density gas. In many cases the cores 
already show central concentration indicative of gravitational 
collapse, and must be associated with the earliest stages of star 
formation. 

We detect 91 dust cores with an 850 fim flux above 
150 mJy/beam (> 4<x detections). These are listed in Table lA. fl 
About half of these are new submm detections. We see from 
Fig. ^ that most of the cores lie in five main groups: from 
east to west in the cloud, the HH211 region, Bl, NGC1333, 
L1455 and L1448. Additionally, there are a few scattered iso- 
lated cores and two cores in B5. These groups are shown in 
more detail in Fig. |2] where the individual cores and the fil- 
amentary structure connecting them are more evident. A full 
inventory of the pre/protostellar population of Perseus includ- 
ing estimates of individual core masses will form the basis of 
a later publication, but it is clear from the submm fluxes that 
these groups contain up to a few solar masses of stars. 

Lower flux sources appear slightly less clustered than 
stronger sources (78% of sources with Fsso^m < 250 mJy 
are in clusters compared to 88% of sources with Fssofim > 
500 mJy), but this may be due to our core selection tech- 
nique, which is less sensitive to weak cores in clusters (see 
Sect.[5)- While there may be additional physical reasons why 
higher mass (and therefore brighter submm) sources preferen- 
tially form in clusters, we are not yet in a position to test this. 

Making a count of cores, we find that more than 80% of 
the cores fall into five main clusters. We define a core to be 
in a cluster if it has 2 or more neighbours within 0.5 pc. With 
our sensitivity to all cores above 0.4 M , this definition guaran- 
tees a mass surface density of these regions of at least 1 .2 M 
within 1 pc 2 . This definition is then consistent with the Spitzer 
tidal stability criterion that a stellar density > 1 .0 M pc~ 3 is 
required i n order for a cl uster not to be disrupted on a timescale 
~ 10 8 yr JSpitzeil[l958l) . although during star formation clus- 
ters almost certainly rely on the gravitational potential of the 
molecular gas for binding. The remaining few isolated cores 
are scattered throughout the molecular cloud. 80% is probably 
a lower limit on the fraction of star formation in clusters in 
Perseus because we underestimate the number of cores in clus- 
ters due to confusion (see below). Infrared studies of embed- 
ded clusters of pre-main-sequence st ars estimate the fr action of 
stars in clusters as 60-9 0% (L1630:lLada et alJll99lh and 50- 
100% (4 nearby clouds: lcarpenterl2000l) . which are consistent 
with our findings. 

This cluste red mode of star formatio n contrasts with the 
Taurus region jAndre & Montmerlel 19941) . where isolated star 
formation dominate s. With the majority of stars in the Galaxy 
forming in clusters jLada & Ladal2003l) . Perseus, like Orion B 
(Mitchell et al]l200ll) is more likely representative of typical 
star formation in our Galaxy. 

3.1. Core clusters and stellar clusters 

Two of the groups of cores are spatially associated with young 
stellar clusters: the HH21 1 group can be associated with IC348, 



and the NGC1333 cores with the NGC1333 cluster, and may 
be evidence for continuing star formation in these clusters. We 
now consider the evolution of star formation in these regions 
by comparing the number of stars in the clusters with the likely 
stellar yield of the current protostellar cores. 

As our pre/protostellar core sample is incomplete below 
0.4 M Q , we take this as the lower mass limit for the count. The 
yield (in stars above 0.4 M Q ) of the protostellar cores can only 
be a rough estimate because of the unknown multiple star frac- 
tion and the unknown final masses of the accreting stars. We 
assume each core yields between 1 and 2 stars above 0.4 M ; 
this takes into account complex cores forming multiple stars. 

In IC348, the infrared population is consistent with con- 
tinuous star formation at a rate of 50 stars with masses above 
0.4 M per Myr from 3.5 Myr to 0.5 Myr ago jMuench et alJ 
I2003L based on an estimate that 2/5 of the total of 348 stars 
shown on the IMF in their Fig. 16, have masses above 0.4 M ), 
This value and a constant star formation rate, suggests ~ 25 
embedded protostars would be expected, if the embedded phase 
has a lifetime of 0.5 Myr. This is consistent with the 18 submm 
cores, containing an estimated 18-35 protostars, we detect in 
the HH211 region. 

Taking an age spread and population for NGC 1333 of 0.5 - 
1 .5 Myr and 143 pre-main-sequence stars JWilking et al.l2004 . 
a constant star formation rate would predict ~ 70 protostars in 
the last 0.5 Myr, whereas we observe 36 cores containing an es- 
timated 36-70 protostars. Here the submm population appears 
consistent or a little less than that expected from a constant star 
formation rate and an embedded phase lifetime of 0.5 Myr. 

An embedded phase lifetime of 0.5 Myr is consistent with 
the youngest T Tauri sources in IC348 and NGC 1333, but 
is longer than estimates for the embedded ph ase lifetime in 
p Ophiuchus and Taurus-A uriga, 0.1-0.4 Myr JWilking et alJ 
Il989t iKenvon et alJll990h . The age spread of T Tauri pop- 
ulations is notor iously difficult to estimate from isochrones 
llHartmannll200ll) . and the youngest IC348 and NGC1333 in- 
frared members may well be less than 0.5 Myr old. However, 
an embedded phase lifetime as short as 0. 1 Myr would require 
either a significant correction to the T Tauri ages or that we 
happen to be looking at a star formation burst after a signifi- 
cant gap in both IC348 and NGC1333, which seems unlikely. 

In conclusion, in IC348 and NGC 1333 our observations are 
consistent with a steady or reduced star formation rate over the 
last 0.5 Myr, but not an increasing one. 

The small population of cores in many of the clusters in 
Perseus raises the question of the cluster IMF - how many stars 
form in large/massive clusters compared to small ones? The 
massive NGC1333 cluster, which has a ste llar population of 
143 stars with a total stellar mass of 79M G JLada et alii 19961) 
contains 40% of our pre/protostellar cores. A further 15% lie 
in the HH211 region bor dering IC348, which has > 300 stel- 
lar members, and 160 M^ JLuhman et alJ2003h . So, in Perseus, 
40-60% of cores are in massive clusters having > 50 M . On 
the other hand small clusters and isolated objects also con- 
tribute 40-60%, much more than the 10% suggested from the 
infrared studies: lLada & Ladal d2003l) suggest that 20-50 M 
clusters contribute no more than 10% of all stars. Unlike our 
statistics on cores, our statistics on clusters are derived from a 
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very small sample of 5, so are less certain. Nevertheless, al- 
though the infrared and submillimetre observations are sen- 
sitive to simila r ranges of masses of stars (> 0.1-0.3M o ; 
ICarpenterkOOoT) . the difference in age between the submillime- 
tre sources and the infrared visible stars may explain this dis- 
crepancy. It is quite likely that the small protostellar clusters 
are not stable and therefore do not classify as clusters by the 
Lada & Lada criteria (> 35 members obeying the Spitzer tidal 
stability criterion) by the time they have become visible at near 
infrared wavelengths. At a typical velocity of 1 km s _1 a star 
moves 1 pc in 1 Myr or 10' at the distance of Perseus. A small 
group of < 10 prestellar cores which evolved to form a 1 Myr 
pre-main-sequence population scattered over a volume 20' in 
diameter could well remain unidentified as a cluster in the in- 
frared. Therefore infrared estimates of the amount of star for- 
mation in small groups may be biassed low. On the other hand, 
in the submm, by taking a snapshot at the earliest stages of star 
formation, we underestimate the count of objects in clusters 
because future star formation is not counted. The smaller clus- 
ters may yet go on to form many more stars. Specifically, Bl 
has a large reservoir of molecular gas and no known population 
later than Class I, and so may be at the earliest stages of cluster 
formation. 

The total mass of gas above a column density threshold of 
6 x 10 22 cm -2 in the clusters and cores is estimated from the 
dust emission to be ~ 2600 M . This assumes isothermal 12 K 
dust: the mass will be lower if the temperature is higher, as 
is the case around the protostars. Lowering the flux threshold 
from the 5<x level to include more low column density mate- 
rial would result in a higher mass. A 5cr cutoff is used to avoid 
contributions from artefacts on large spatial scales due to the 
reconstruction of the chopped map. With the 5cr cutoff, this 
contribution is limited to be less than 30% of the total: the mass 
estimate for the clusters alone (NGC1333, HH21 1, Bl, L1448, 
L1455) is 1800 M Q . The mass fraction of the molecular cloud 
actively involved in star formation at this time is less than 20% 
of the tota l gas mass of 17,000 Mp estim ated from the visual 
extinction dBachiller & Cernicharoll986bl) . The mass of dense 
molecular gas traced by C 18 1-0, which has a column density 
threshold for detection of 1.6 x 10 21 cirT 2 (A v =* 2) and a crit- 
ical density of 400 cirT 3 , is ~ 6000 M Q , lower than the mass 
estimate from extinction as it excludes low density and atomic 
gas. Thus a greater fraction of high column density gas than 
low column density gas is in the active star formation regions, 
up to as much as half of the gas above A v 2. 

4. Filaments 

The SCUBA 850 /mi image shows filamentary structure 
throughout the cloud (Fig. 0. The active protostellar clusters 
are linked to dense filaments: the horseshoes of HH211 and 
Bl, the arcs of L1455 and L1448, and the zigzag of dense gas 
in NGC1333. Weaker filaments imaged by SCUBA to the east 
of HH211, the south of Bl, the south of NGC1333 and around 
L1448 each reach up to 1 pc in length, with the upper limit here 
presumably set by our sensitivity limit. The dust filaments have 
typical full width half maxima of 1'. Filaments to the east of the 
HH21 1 cluster and to the southwest of Bl do not appear to be 



connected to current protostellar activity but may be fragment- 
ing into cores to form the next round of stars. All the filaments 
show intensity fluctuations along their length, with the most 
fragmented filaments already identified as strings of cores. 

On much larger scales, ~ 20 pc, the C ls O map (Fig. 
which traces material down to lower column densities, also 
shows filamentary structure. The dust filaments detected by 
SCUBA are found in regions of high C ls O column density. 
In fact the entire Perseus cloud can be described as con- 
sisting of one fi lament with a continuous velo city structure 
along the cloud dBachiller & Cernicharolll986al) . SCUBA de- 
tects 850 fim emission at the 5-<x level at H2 column densities 
of ~ 6 x 10 22 crrT 2 (A v - 60), as the C ls O starts to become 
optically thick. This column density is rarely reached and the 
SCUBA map is mostly empty. Some extended structure is also 
lost in the reconstruction of the SCUBA maps because of the 
limited range of chop throws. 

The filaments a re strongly reminiscent of str uctures formed 
in turbulent flows dMac Low & Klessenll2004l) . What we see 
in the dust continuum is that the filaments are observed down 
to the small scales and high densities associated with cluster 
formation. This suggests that not only do molecular clouds on 
the parsec scales of the whole Perseus molecular cloud form 
by turbulence, but that the role of the filaments persists down 
to the cluster formation scale. A more detailed investigation of 
the role of turbulence, including kinematical information from 
molecular lines, is left to a future paper. 

The mass in the filaments is substantial with a mean mass 
per unit length of between 4.7 and 1 1.5 M per 0.1 pc (~ 1'), 
assuming the filaments are isothermal at 12 K, and taking a 
sample of 43 filament cross-sections. The uncertainty is due to 
the subtraction of background flux, some of which is extended 
flux belonging the filaments and some of which is an artefact 
of the reconstruction process. An estimate of the artificial back- 
ground in regions of blank sky suggests that the true filament 
mass lies roughly midway between these values. 

The filament mass per unit length of 5-11 M e per 0.1 pc 
(1') is 2-4 times the traditional thermal Jeans mass for a Jeans 
length of 0. 1 pc but this does not imply that the filaments have 
to be unstable. The cr iteria for a fila ment are different for those 
of a uniform medium dLarsonll985t) . Theoretical models for fil- 
aments range fr om a simple non-magnetic isothermal filament 
dOstrikerl 1 1964b to magnetised filaments w ith helical fields 
which can be toroidal or poloidal or both dTillev & Pudritzl 
l2003UFiege & Pudritzll2000HFiege et al.ll2004 . A kev predic- 
tion from these models is the ratio of the mass per unit length 
m to the cylindrical virial mass per unit length m wn , which is 
given by 



where cr tot is the 1 -dimensional total (thermal plus non- 
thermal) velocity dispersion of the average gas molecule. For 
filaments with or without magnetic field, 

Wvir\ \W\) (Py 

where M/\W\ is the ratio of ma gnetic energy to gravitat ional en- 
ergy per unit length. Following lFiege & Pudritzl ( 120001) Table 2, 



8 



Hatchell et al.: SCUBA Perseus survey 



we assume the external pressure to be 10 4 5 K cm 3 , and esti- 
mate the internal pressure from 



P/k 



2x 10 6 KcirT 



using an average density in the filaments of 10 5 cirT 3 . The exact 
values of these quantities are unimportant; what is important is 
that the ratio of external to internal pressure P$/(P) — 0.01, 
which is much less than 1, so there are three regimes: m/m vlI ~ 
1, the filament can exist without magnetic support, m/m vlI » 1 
a supporting poloidal field is required, and m/m vn <s 1, when 
a binding toroida l field is required (see lFiege & Pudrital2000t 
iFiege et al.l2004l for details). 

From a sample of 7 C 18 1-0 spectra towards filaments, we 
measure velocity FWHM of 0.6-1.0 km s _1 . The conversion 
from velocity FWHM Av b s measured in C ls O to <T tot , taking 
into account that the C ls O molecule is more massive than the 
predominant H2, is 



2 



1 / , kT ( 1 1 

AvL + 81n2— -- 



81n2 



mp\m m C i8o 



where nip is the proton mass, m — 2.33 assuming 1 He for every 
5 H2 molecules, T — 10 K, and mciso = 30. 

The measured C ls O linewidths correspond to virial masses 
of 45-100 M e pc~', compared to the measured mass per unit 
length of 47-115 M pc~'. This simple analysis suggests that 
the Perseus star formin g filaments are consistent with being 
non-magnetic filaments ( lOstrikeril9 64). 



This last result disagrees with the 



Fiege & Pudritzl feOOd) 
studies of star-forming filaments, which typically used 13 CO 
linewidths, and found mjm^ lx < 1. From a sample of 8 13 CO 1— 
spectra, we measure velocity FWHM of 1.4-2.0 km s _1 , typ- 
ically a factor of 2 greater than C ls O linewidths. If instead 
these 13 CO linewidths are representative of the velocity field 
within the filaments then ours too require additional binding by 
a toroidal field. However, it seems more likely to us that the 
C I8 linewidths are more appropriate for the thermal support 
within the dense filaments, as C 18 profiles are less affected 
by optical depth, and that the correct solution is one where the 
total toroidal field is small. 

The mass derived from the dust emission, which increases 
if either the dust temperature or the 850 /mi emissivity is lower 
than those assumed, both of which are possibilities, adds an- 
other uncertainty. An examination of the transverse column 
density structure and the dust polarisation would provide a fur- 
ther test of the filament support but this is beyond the scope of 
the current paper. 

The alternative conclusion that could be reached from the 
high virial mass is that the filaments are not stable but tran- 
sient, dispersing on the crossing time of 0.1 Myr (assuming 
a filament width of 0.1 pc and velocity dispersion from C 18 
of cr tot = 0.4 km s _1 ). If this is true, then the star formation 
which is clearly occurring in the filaments must also happen on 
a crossing time. 

Given filament masses of 5-1 1 M per 0. 1 pc, there is no re- 
quirement for substantial large scale flows along or onto the fil- 
aments in order to gather sufficient material for further star for- 
mation. Typically, in regions where cores have already formed, 



one or two cores are found within 0.1 pc, and local fragmen- 
tation will suffice to explain their masses. On the other hand, 
each filament does not hold a vast reservoir of mass - a few 
tens of M Q . Therefore the filaments are not required to dis- 
perse large fractions of their mass in order to remain consistent 
with the observed stellar densities in IC348 or NGC1333. In 
other words the conversion of dense molecular gas into stars 
may be fairly efficient once the filament stage is reached, and 
most of the filament mass may ultimately go into cores. This 
potentially high conversion efficiency of dens e molecular gas 
into s tars has also been noticed for NH3 cores jFuller & Mversl 
119871) . Alternatively, stellar clusters may be built up by a se- 
ries of filaments that form a few stars inefficiently, disperse, 
and are replaced, requiring ongoing filament production in the 
same location. 

5. Column density condition for star formation 

We are interested in the conditions under which stars form and 
it is clear that the density of the parent molecular cloud is an 
important factor. Density itself is not so easy to measure di- 
rectly from dust continuum or C 18 emission, but what we can 
consider is how the numbers of stars formed varies with the 
line-of-sight column density (or equivalently visual extinction). 
An extinction threshold for star forma tion at A v - 7-9 was 
found in both Tau rus and Ophiuchus Jjohnstone et alJ Eo04l: 
lOnishiet all 19981) . which could indicate that a certain amount 
of shielding from the i nterstellar rad i ation f ield is required for 
star formation to occur. lOnishi et al] l l 19981) measured the col- 
umn density of C ls O associated with sources in Taurus, and 
found that the cold IRAS sources (assumed to be the younger 
protostars) and starless H 13 CO + cores had a average column 
density of at least 8 x 10 21 cm 2 or A v = 9. Ijohnstone et alJ 



( 2004) compared 850 yum cores in p Ophiuchus with an opti- 
cal/IR extinction map and found no submm cores below an A v 
of 7. In both cases the interpretation is that star formation is 
inhibited below a certain A v . 

The C 18 integrated intensity is a good tracer of total gas 
column density over the typical range of column densities in 
molecular clouds, but is less reliable at extreme column den- 
sities. At low column densities the C 18 abundance falls due 
to photodissociation. On the other hand at high densities the 
line becomes optically thick and in addition in low tempera- 
ture, hig h density regio ns the C ls O depletes onto dust grains 
(e.g. lTafalla et all2004 and so ceases to trace the gas. But be- 
tween these e xtremes, a linear C 1 8 Q - A v relationship holds 
above A v ~ 4 dFrerkingetalJl982l) . 

We see in Fig.^that the C 18 is relatively smooth on scales 
of a few arcminutes across the cloud, rather than peaking at the 
submm cores. C 18 integrated intensity at 1' resolution is a 
poor tracer of submm cores, tracing instead the environment 
on larger scales. This is probably largely due to beam dilution 
- the C 18 effective beam area is (l'/14") 2 = 18 times greater 
than the 850 fim beam, so a core of as much as 0.7 M Q could 
produce less than 1 K km s in C I8 integrated intensity if it 
had no extended envelope. Depletion also helps reduce the con- 
tribution of core emission to the C I8 0. C ls O deple tion has been 
measured in prestellar cores on scales of 0.1 pc JTafalla et alJ 
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Table 1. Fraction of cloud by area and number of cores in each 
C ls O integrated intensity range. 
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0.2 
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15.9 


> 5 


> 19.5 


0.04 


1 


1.1 



1 Fraction of the area of the cloud mapped in C 18 



12004 . with the C l8 abundance substantially reduced in the 
dense parts of the cores. Finally, the C 18 1-0 line saturates at 
about the same column density at which the SCUBA 850 /mi 
emission becomes visible, which further reduces the contribu- 
tion from the highest column density lines-of-sight. The result 
is that the C 18 traces the environment of the cores, and thus, 
if we look at young enough cores, the environment in which 
cores form. 

5.1. Probability of finding a core 

In this section we investigate how the probability of finding a 
core varies with the C 18 integrated line intensity on arcminute 
scales. 

Table ^ compares the area of the cloud in each C 18 in- 
tegrated intensity range with the number of submm cores in 
that region. Although there are clearly cores associated with 
areas of the cloud with a wide range of cloud column densi- 
ties, cores are much more likely to be found in regions of high 
column density. In Fig. [5] we demonstrate this graphically by 
plotting the probability of finding a submm core, P(core), as 
a function of integrated C ls O intensity, 7(C 18 0) (K km s -1 ). 
The probability of finding a core is derived from two observed 
distributions. The first is the C ls O integrated intensity towards 
the peak of each submm core, determined from the C ls O map 
(Fig. \^p). We bin the results by C 18 integrated intensity to 
generate N COKS (Ii), the number of submm cores in each inte- 
grated intensity bin. The second is a count of the number of 
pixels in the map in each integrated intensity bin, A^ p i xe is(/i). 
This gives a measure of the (probability) distribution of inte- 
grated intensities P(7) in the surveyed cloud, or in other words 
how much of the cloud area lies in a particular integrated 
intensity range. The probability of finding a submm core at 
a position with C ls O integrated intensity I is then given by 

P(COre) = jVcoresW/NpixelsW- 

As both distributions N cores {Ij) and N p i xe \ s (Ii) at some point 
contain very small numbers of counts we have to be very care- 
ful with the error budget in calculating this quotient. We take 
the Bayesian approach detailed in Appendix |A| to calculate 
P(core). 

Fig-E]shows that P(core) is a steeply increasing function of 
integrated intensity I. In other words, the probability of find- 
ing a core increases rapidly with C I8 column density. Fitting 




|"t ub (C 1b 0) dv / K km s" 1 



Fig. 3. Probability of finding a submm core in a 22" square 
pixel with a given 1' beam-averaged C 18 integrated inten- 
sity: observed (thin black line)and probability distribution of 
log 10 (P(core)) = N C0TSS (Ii)/Np ixels (Ii) for each integrated in- 
tensity bin (greyscale), and best-fit power law (thick black 
line). To convert the probabilities to another pixel size use 

n, 7i/ x (pixel area) 
P' = P(core) x — — 7pp_ -. 

(for simplicity, and for want of a theoretical prediction) a power 
law kl° to the distribution gives a best fit result of a — 3.00^ 8 , 
logiok = -3.15+0.08 (uncertainties are 68.3% confidence lim- 
its), k is the scaling factor which gives the probability of finding 
a core in a 22" square pixel. A power law must break down as 
a description of P(core) at high 7(C 18 0) as the probability of 
finding a core in any one pixel approaches 1 but this limit is not 
approached with our pixel size. 

The exact shape of the distribution is uncertain at the low 
end because of incomplete map coverage and lack of sensitiv- 
ity. Both JCMT and FCRAO maps have limited coverage of 
low column density regions, and as the area mapped in C ls O is 
larger, we may have missed SCUBA cores. So we may underes- 
timate (by only a few cores, we believe) the probability of find- 
ing a core below 7(C 18 0) = 1.0 Kkms -1 . Fig. □ also shows 
that the C I8 observations were not biassed towards regions 
with SCUBA clumps, so we do not overestimate the number 
of cores at low intensity. The measurement uncertainties also 
become significant below 7(C ls O) = 1 Kkm s -1 (A v = 5.2) - 
the C I8 RMS is 0.3 K km s _1 . Therefore, the source count in 
the lower 5 integrated intensity bins (Fig. O is also uncertain. 
We made a fit excluding the lowest column density channels 
(7(C 18 0) < l.OKkms -1 ), which steepens the power law to 
a = 3.45 + 0.27, a change which is significant only at the 95% 
level. 

What does P(core) oc 7 3 mean? We can rule out an abso- 
lute threshold in N( C 18 0) (and by implication A v ) for submm 
cores (though not necessarily for star formation). There are 
clearly real submm cores at low C ls O integrated intensities of 
4 K km s or visual extinctions as low as A v = 3 (N(H2) 
2.7 x 10 21 cm -2 ). At these column densities the uncertainty on 
the C ls O data is high but the low extinctions are independently 
confirmed by 13 CO data from the same observing run. There 
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are submm cores at very low C O column density despite the 
fact that the peak column densities in the cores must be high 
to be detected with SCUBA, and that the observations are less 
sensitive to the larger, more diffuse cores which form in low 
pressure regions (Sect.|3Jl. 

The 10% of cores with I < 1.0 K km s _I include 
the L1455 cluster plus IRAS 03235+3004, 03422+3256, and 
03262+3123. The question remains whether stars can form 
at these low column densities or if these low column density 
cores are all more evolved (Class I) protostars, or whether the 
lower angular resolution of the line data is biasing the column 
density estimate to a low value. Certainly the majority of the 
/ < 1 K km s _1 (A v < 5.2) sources have IRAS identifications 
and are therefore likely to be classified as Class I. Potentially 
less evolved cores are L1455 smml (ff(2000) = 03 h 27 m 42 s :9, 
5(2000) = +30°12'28"), which has a CO outflow, is definitely 
protostellar, and has low IRAS fluxes, L1455 smm2 (a(2000) = 
03 h 27 m 46 s i6, 5(2000) = +30°12'05"; probably starless), and 
three weak cores to the w est of B 1 which may be noise arte- 
facts. |Ojnshi_et al. ( 1998) found that cold cores in Taurus were 
only found at a high C 18 threshold of 8 x 10 21 cm 2 , equiva- 
lent to /( C ls O) ~ 2.0 K km s _I or A v = 9, whereas warm cores 
(more evolved IRAS sources) could appear at lower C ls O col- 
umn density. This appears to roughly hold in Perseus as well 
as Taurus, but we need to investigate this further once we can 
better classify the Pers eus sources. All the s ubmm cores in 
Ophiuchus lie at A v > 7 Jjohnstone et al.l2004l) . As the cores in 
Ophiuchus were identified in the same way as in this study of 
Perseus, that is, from SCUBA 850 fim observations, there must 
either be a real difference between star formation in the two 
clouds, or a significant difference in the way that optical/IR ex- 
tinction and C 18 trace column density in the range A v = 2-7. 
p Oph is lacking in starless cores and Class protostars com- 
pared to other dark clouds and there may have been a burst of 
star formation 10 s yea rs ago in this region that has now tailed 
off JVisser et al.ll2002t) : in Perseus, however, we have clear ev- 
idence for continuing star formation, so there are differences 
between the ages of the populations in the two clouds which 
might underly the differing column density thresholds, though 
the mechanism is unclear. 

The 7 3 power law rules out simple models in which the 
number of cores is proportional to the mass along the line- 
of-sight (P(core) <x I). If density is the determining quantity 
for star formation, central densities must rise steeply with in- 
creasing column density. Theoretical models which reproduce 
the observed C 18 integrated intensity - core probability rela- 
tionship are required to find the actual physics behind the 1 3 
power law (if indeed a power law is the right functional form). 

This steep rise in P(core) with column density is consistent 
with the overall trend seen in star forming regions for high col- 
umn density regions to be associated with formation of larger 
numbers of stars. Isolated low mass star formati on in Taurus 
occur s at column densities of = 8-15 x 10 21 cirT 2 JOnishi et alJ 
Il998l) . Column densities for clustered star formation in Perseus 
reach 2 x 10 22 cirT 2 . In contrast, column densi ties in massive 
star fo rming regions typically exceed 10 23 crrT 2 JHatchell et alJ 
119981) . 



We have already mentioned C ls O depletion and optical 
depth, both of which act to reduce the observed /( C 18 0) for a 
given H2 column density. If plotted against H2 column density, 
the P(core) distribution would flatten at high column densities. 
A rough correction for the optical depth suggests that the power 
law index could reduce from 3.0 to ~ 2.0 in the worst case. This 
is still significantly steeper than linear. 

Given a C 18 map and the above P(core) vs. /( C ls O) rela- 
tionship, we can now predict - at least statistically - how many 
submm cores will be found and how they will be distributed. 
It will be interesting to see how well this relationship holds in 
regions other than Perseus. There are obviously more factors 
involved in core formation than C 18 column density, and a 
better model of where stars form will need to take into account 
filaments, clustering and other factors. 

6. Conclusions 

We have mapped the submm dust emission from the Perseus 
molecular cloud with SCUBA at JCMT with 5000 AU resolu- 
tion, revealing 91 embedded starless and protostellar cores and 
the filamentary structures which lead to star formation. The to- 
tal mass of the Perseu s cloud derived from extinction measure- 
ments is 17,000 M Q JBachiUer & Cernicharolfl986bh . Of this, 
6000 M is traced by C 18 and 2600 M e by dust in cores and 
clusters. We conclude: 

1. Only a small fraction, less than 20%, of the mass of the 
molecular cloud is involved in star formation at this time. 

2. By number, 80% of the star formation is occuring in clus- 
ters, defined by a stellar density of more than 1 M pc -3 . 
Some current star formation activity may be associated 
with the larger stellar clusters IC348 and NGC1333, but 
about half the star formation is taking place in small clus- 
ters containing less than 35 members. 

3. Filamentary structure is evident both in the molecular gas 
on 20 pc scales and on small scales in the dust emission. 
High column density dust filaments with masses per unit 
length of 5-11 M Q per 0.1 pc are associated with the ex- 
isting clusters. This range of masses is in agreement with 
models of filaments with either u nmagnetised filaments or 
a slightly dominant toroidal field (Fiege & PudritzbOOol 

4. The probability of finding a submm core is a steeply ris- 
ing function of column density (as measured by C 18 inte- 
grated intensity): P(core) oc 7 3 . There is no C ls O column 
density cutoff (and by implication A v ) below which there 
are no cores. 

These results are specific to Perseus, and it will require 
similar studies of a much larger sample of such clouds before 
we can fully understand which of these features are typical of 
star forming clouds and which are, for some reason, unique to 
Perseus. 
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Appendix A: Details of power law fitting procedure 

In this Appendix we give details of the analysis of the core probability 
vs. C ls O column density. We first assume that the counts in each inte- 
grated intensity bin /,, x, = iV cores (/,) and y, = iVpjxeuf/;), obey Poisson 
statistics, i.e., are related to the true values fi x and ji y by: 



- f(y) 



v(y\th) 



/4 e 



By Bayes' theorem we can reverse this. We take a a flat prior with 
a positivity constraint on fi x and /j T , PQu A -) = 1 (ji x > 0) or (pi x < 
0) and similarly for [X y (see, e.g., ISiviall 19961 for an introduction to 
Bayesian data analysis). Then 



PQi y \y) =P(y| Mv )xP(^): 



l4 e 



(A.l) 
(A.2) 



and similarly for fly. 

We calculate the data-derived probability of finding a core at inte- 
grated intensity /,, C b s ,/ = x,/y, directly. However, this is only an esti- 
mate of the true value of C because of the uncertainties in x, and y,-. We 
next calculate the PDF for the true value of finding a core P(C Ul i e ] x, y) 
given Xj and y f . To obtain the probability distribution of the true quo- 
tient Ctruc given the observed x and y we rewrite the x distribution 
(Equation lA, U in terms of fi y and C lnre : 



P(yU Y | X, C^e) 



(Ctme)^ 



Cuucfly 



(A3) 



Now with two probability distributions on fi y (Equations IA.2I 
and lA.3> we can multiply and integrate out fi y : 



P(C lmt \x,y) oc 



Jo y! 



dfi y 



x (x + y)^ Cj ae 
(1 + C trac ) (j ^ +1) x\y\ ' 

Finally to find the most probable power law fit, for each /, we 
calculate P(C tmc , i = kl" | x ; ,y;) from the previous formula and 



P(k,a\x,y) = Y\v(k,a\ x b y t ). 



The PDFs for C lllle are plotted as greyscales and contours for each in- 
tegrated intensity bin on Fig. [3] We marginalise (integrate over) k and 
a in turn to estimate the uncertainties in each parameter individually 
as given in Sect.|5] The PDF for a is shown in Fig. lA. II 
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Table A. 1 : SCUBA 850 //m core positions, names and references, peak 850 yum 
fluxes and C ls O integrated intensity at the submm peak position. 



No. 


RA2000 
hh:mm:ss 


D ec 2000 

dd:mm:ss 


Cluster* 1 1 


Name 


Reference 1 ' 


FTn 

1 850 

mjy/bm 


/(C I8 0) (4) 
K km s 1 


1 


03:33:17.9 


31:09:33 


Bl 


Bl-c 


Matthews & Wilson (2002) 


1638 


4.0 


2 


03:33:21.4 


31:07:31 


Bl 


Bl-bS 


Hirano et al. (1999): Matthews & Wilson (2002) 


1554 


4.8 


3 


03:33:21.1 


31:07:40 


Bl 


Bl-bN 


Hirano et al. (1999): Matthews & Wilson (2002) 


1200 


4.8 


4 


03:33:16.3 


31:06:54 


Bl 


Bl-d 


Matthews & Wilson (2002) 


514 


4.1 


5 


03:33:01.9 


31:04:23 


Bl 


None 


218 


3.3 


6 


03:33:05.6 


31:05:06 


B 1 


None 




207 


2.8 


7 


03:33:16.5 


31:07:51 


Bl 


Bl-a 


Hirano et al. (1999): Matthews & Wilson (2002) 


202 


4.5 


8 


03:33:04.0 


31:04:56 


Bl 


None 


202 


2.9 


9 


03:33:03.0 


31:04:38 


Bl 


None 




196 


3.4 


10 


03:33:27.3 


31:07:10 


Bl 


None 




179 


4.2 


1 1 


03:33:00.3 


31:04:17 


Bl 


None 




160 


3.4 


12 


03:43:56.5 


32:00:50 


HH211 


HH211 


McCaushrean et al. (1994): Chandler & Richer (2000) 


1494 


2.5 


13 


03:43:56.9 


32:03:05 


HH211 


IC348 MMS 


Eisloffel et al. (2003) 


1032 


3.1 


14 


03:44:43.9 


32:01:32 


HH211 


None 


623 


2.0 


15 


03:30:15.5 


30:23:43 


HH211 


No SMM/MM source 




358 


2.2 


16 


03:44:01.0 


32:01:55 


HH211 


None 




296 


2.4 


17 


03:43:57.9 


32:04:01 


HH211 


None 




247 


3.0 


18 


03:44:03.0 


32:02:24 


HH211 


No SMM/MM source 




237 


2.5 


19 


03:44:36.8 


31:58:49 


HH211 


None 




234 


1.8 


20 


03:44:05.5 


32:01:57 


HH211 


No SMM/MM source 




209 


2.4 


21 


03:44:02.3 


32:02:48 


HH211 


None 




207 


2.9 


22 


03:44:06.2 


32:02:12 


HH211 


None 




204 


2.4 


23 


03:43:37.8 


32:03:06 


HH211 


None 




170 


1.5 


24 


03:43:42.5 


32:03:23 


HH211 


None 




168 


2.3 


25 


03:44:48.5 


32:00:31 


HH211 


None 




153 


2.0 


26 


03:43:44.4 


32:02:56 


HH211 


None 




153 


2.5 


27 


03:25:35.9 


30:45:30 


L1448 


L1448NW 


Loonev et al. (2000): Barsonv et al. (1998) 


2300 


2.9 


28 


03:25:36.4 


30:45:15 


L1448 


L1448N B 


Loonev et al. (2000): Barsonv et al. (1998) 


2243 


2.7 


29 


03:25:38.8 


30:44:04 


L1448 


L1448 C 


Bachiller et al. (1991): Barsonv et al. (1998) 


1737 


2.3 


30 


03:25:22.4 


30:45:11 


L1448 


L1448-IRS2 


O'Linser et al. (1999) 


1316 


2.2 


31 


03:25:25.9 


30:45:03 


L1448 


None 


429 


2.4 


32 


03:25:49.0 


30:42:25 


L1448 


None 




364 


2.0 


33 


03:25:50.6 


30:42:01 


L1448 


None 




306 


2.1 


34 


03:25:30.8 


30:45:07 


L1448 


None 




231 


2.6 


5 J 


Uj.z/.jy.l 


3U. 13.U1 


T 1 A £^ 


KJNU 1j rlK 


Kensel et al. (zUUz) 


JJO 


1 1 

1 .4 


36 


03:27:42.9 


30:12:28 


L1455 


None 


402 


0.9 


37 


03:27:48.4 


30:12:09 


L1455 


No SMM/MM source 




324 


0.8 


38 


03:27:46.6 


30:12:05 


L1455 


No SMM/MM source 




272 


0.9 


39 


03:27:38.1 


30:13:57 


L1455 


No SMM/MM source 




238 


1.4 


40 


03:27:39.9 


30:12:10 


L1455 


None 




232 


1.4 


41 


03:29:10.4 


31:13:30 


NGC1333 


IRAS 4A 


Sandelletal. (1991) 


7000 


3.7 



Table A.l: continued. 



INO. 


P A 

hh;mm!ss 


J-*cC2000 

dd:rnm:ss 


L-iuster 


Name 


Reference^' 


^850 

mjy/bm 


r/v 1 

^(,*-18^) 
IZ 

JV Kill S 


A O 

42 


07 .70. 1 7 o 

Uo.2y. 12.U 


7 1 . 1 7 . 1 
31.13. 1U 


JNUL.IOOO 


IO AC /ID. C \Z 1 

lKAc> 4d, ioJv-o 


c .. ,-,^1 1 at „i /iooi\. c . . , 1 ,^ 1 1 p t iz ,, /onoiK\ 
ioanaell et al. ( IWI ), oanaell <sc Jvnee (2UU1D) 


707 C 

oU2o 


7 /I 
0.4 


■40 




7 1 ■ 1 S- C Q 

1 . 10. J7 


wnr 1 1777 

1 N V J v 1000 


11117 11 \/T\/TQ 10/17 
nn /-ii iviivio i,oiv-iz/io 


V_I11IU Cl al. \ 1 77 / ) 


087 1 
ZOO 1 




1 1 

44 


Uo.2o. jj. o 


71-1/1 . 7A 
Ol. 14.00 


Mf'Y" 17 77 
JNOL-lOOO 


1 AC OA. C 1Z O 

IKAo 2A, 0K-0 


Onn/lall at „1 /'^Ci(\A\• C ., A ,^ 1 1 O r L^naa /"OOOIKA 

oanaell et al. ( 1W4), oanaell & Jvnee (2UU1D) 


77 C C 


2 
0.0 


4j 


uo.zy.ui.4 


7 1 -70.70 

0i.2u.2y 


M^f 1 777 
INLjC 1 000 


CT7 0/1 

olv-24 


^anueii oz Jvnee (zuuioj 


I 1 1 o 

I I iy 


7 7 

0. / 


1 A 

40 


07 .70. 1 o n 

uo\2y. iu.y 


7 1-1 .77 

1. lo.Z / 


MfP 1777 

JNUL-looo 


oK-2U 


ioanuell oz Jvnee (2UU1DJ 


i nA 1 
lUol 


7 A 
0.0 


4 / 


ni.io.<n 7 
Uo.zo. jV. / 


7 1 -71 .7/1 

ol. 21.04 


M/" 1 /" 1 17 77 

JNOL-looo 


olv-o 1 


ioanueJl oz Jvnee (2UUJD) 


770 

/ /u 


1 o 

4.v 


4o 


07 .70. 1 7 A 

Uo.2y. 1 o.o 


1. lo.oo 


vt/" 1 rim 
JNOL-looo 


TD AC /IT. Cf c 

IKAo 4L., oJv-o 


t ^^.«ai7 at ni /'onnnv c<i«^iaii p T Tz„ a£ , /onmu\ 
.Looney et aJ. (2UUU), banaeJJ & Jvnee (2UUJD) 


AA7 
00 / 


2 2 
J.J 


4V 


O7-70-7A 7 

Uo.2o.oo. / 


7 1-17 .70 

ol. I0.0U 


Mf'Y" 17 77 

JNOL-looo 




Sandell & Knee (2001b) 


AO 1 

02 J 


O 

2.U 


oU 


07.70-OA C 

uo.2y.uo.o 


t 1 . 1 e .on 

ol. I0.0V 


VT/" 1 /" 1 1777 

JNOL-looo 


Till n 1 1 TVTIV/TC /I.CF 1 C 

HH /-ll MJVlo 4,oJv-lo 


/~M".J«J at ^1 /OOm \. Cnn/lall JP T T^"«aa nAf!1K\ 

L.mni et aJ. (2UUJ ), oanueJJ oz Jvnee (2UUJD) 


A 1 Q 

oiy 


1 

0.2 


C 1 

J J 


O7-0Q.no 

uo.2y.uo.o 


li.ic.io 

1. lo. lo 


Mff 1 777 
IN LjL.1 000 


O K- 1 


^anueii oz jvnee (2UUioj 


/4 


7 Q 
2.5 


c o 
02 


07-70.07 7 

uo\2y.uo\ / 


7 1 . 1 /I . C 7 

ol. 14.00 


M/" 1 /" 1 1777 
JNUL.I0OO 


I 1 1 1-7 1 1 A/TA/rC A- C 1 A 

HH/-11 IV1 IV1 0, ioK-14 


Chini et al. (2001): Sandell & Knee (2001M 


C 1 7 

1 / 


7 

o.o 


c ~2 
3 o 


07-70.0/1 C 

uo\2y.u4.o 


7 1 .70. ?n 
0l.2U.oV 


M/T^ 1 1777 
JNUL.I0OO 


oK-20 


Cn«^iaii i? r v«aa /'onmi-.\ 
ioanueJl oZ Jvnee (2UUJDJ 


C 1 7 

o lo 


1 

4.0 


c i 
04 


07 .70. 1 O 7 

Uo.2y. 1U. / 


7 1 -71 . A C 

ol.21.4o 


JNLjL.1ooo 




Onn/lall S? r Vnaa /"OnniK\ 

oanueJJ oZ Jvnee (2UUJD) 


A OO 

42U 


7 A 
O.O 


cc 

J J 


07-7C/10 /I 

Uo.zo.4U. 4 


71 .n.c 1 
ol. 1 /.ol 


JNOL-looo 


None 


/I 1 o 

41U 


1 o 

i .y 


nA 


07-70.07 7 

Uo.zy.u/. / 


7 1-71. C7 

1.21.0 / 


UPf 17 77 

JNOL-looo 


c \z on 
oK-2y 


oanueJJ oZ Jvnee (2UUJD) 


1QA 

oV4 


1 7 

4.0 


J / 


07-7Q- 1 7 

uo.zy. lo.z 


7 1 . < - 1 1 

1 .20. 1 1 


Mff l 777 
InLjC 1 000 


CT7 QO 

olv-00 


^anueii oz jvnee (zuuidj 


77A 
/O 


4. 1 


eg 

JO 


07-70-7/1 O 

Uo.2y.24.U 


7 1 -77.71 
1.O0.21 


M/" 1 /" 1 1777 

JNOL.1 000 


None 


7 7 1 

O / 1 


n o 

u.y 


<o 
oV 


07 .70. 1 A C 

uo\2y. io. j 


7 1 . 1 7 . 7 C 
Ol. 12.O0 


M/" 1 /" 1 17 77 
JNUL.I0OO 


None 




7 7 1 
/ 1 


7 

0.2 


An 
OU 


07-7Q-7Q 7 
Uj.ZO. J7.J 


71-1 0-77 
1 . 10.2 / 


XTP/^ 1 777 

looo 


None 




oo4 


o o 

2.v 


1 


07 .70. 17 7 


7 1 -77. CO 

ol.2 /.oU 


VT/" 1 /" 1 17 77 
JNOL-lOOO 


None 




7/10 

o4y 


1 

2.0 


o2 


07-70.07 1 

Uo.zy.u/. i 


7 1-17.7/1 
ol. 1 /.24 


VT/" 1 /" 1 1777 

JNOL-looo 


C IX 1 

0K-I0 


oanueJJ oZ Jvnee (2UUJD) 


7 10 

o lo 


1 A 

4.U 


A~5 
00 


07-7Q- 1 Q 

Uj.Zy. lo.o 


7 1 - 77- 1 7 
1.20. 1 / 


MPf 1 777 

INOC 1 000 


CT?" 70 


^anueii oz ivnee (zuuioj 


7 no 
oUo 


7 Q 
O.O 


A 1 

04 


07.70. oc c 

U7.zy.zj. j 


7 1-70.10 

ol.2o.lo 


M/" 1 /" 1 1777 

JNLjL.1ooo 


None 




707 

2V2 


1 1 

1 .4 


AC 

Dj 


07.70.no A 
Uo.2y.UU.4 


7 1 - 1 7 . n 1 
1. 12. Ul 


VT/" 1 ^ 17 77 

JNLjL.1 000 




Sandell & Knee f2001M 


700 

2oo 


7 

2.0 


AA 


07-70-OC 7 

Uo.2y.Uo. 


7 1 .77 . 1 1 
01.22. 1 1 


KT/" 1 /" 1 17 77 

JNOL-looo 


C IZ or\ 


loanaeJJ oZ Jvnee (2UUJD) 


770 

22o 


C 1 

J. 1 


A"7 
/ 


07 .70. 1 n 7 
Uo.2y. IV. / 


7 1 -77 . CA 

ol.2o.oo 


Mr/^ 17 77 

JNOL-looo 


None 


777 

22o 


1 7 
4. / 


Do 


Uo:2o:oo.2 


7 1 . 1 n. 1 7 

ol: lv: lo 


JNOClooo 


No SMM/MM source 




OO 1 

2U1 


2 O 

o.y 


AO 

oy 


07. 7C 7/1 /I 

Uo.2o.o4.4 


7 1 .nA. co 
1.U0.0V 


KT/" 1 ^ 1 7 7 7 
JNUL.I0OO 


None 




7 n i 
2U J 


1 7 
1 .0 


7n 
/U 


07 .70. 1 c 7 
Uo.2y. 1 j. 


7 1 .on. 7 1 
0l.2U.oi 


M/^r 17 77 

JNLrL.looo 


olv-22 


ioanuelJ oZ Jvnee (2UUJbj 


1 Q7 

JV2 


1 

4.2 


7 1 
/ 1 


U0.Z0.7O. / 


7 1 -n^-^7 
l.Uo.o / 


Wf^f^ 1777 
InOCIooo 


None 


1 77 

J / / 


1 

i .v 


"7 7 


07 .70. 1 n 1 
Uo.2V. iy. 1 


71-1 1 ■ 70 

ol. 1 1.00 


MPf 17 77 

JNOL-looo 


No SMM/MM source 




1 A7 
10 / 


A 

2.0 


~! 2 

1 J 


07.70.70 
Uo.zo.00.0 


7 1 . 1 n. 1 c 
1. IV. 1 


\ YY" 17 77 

JNOL-looo 


None 




1 AO 

lo2 


1 

2.4 


7 1 

/4 


07.70.77 c 
Uo.zo.oz. J 


71.11 .no 
1. 1 1 .Uo 


Mrr 1777 
JNLjL.1ooo 


None 




1 <^A 

loo 


1 o 

i .y 


~7C 


O7-70-/17 A 

Uo.2o.42.o 


7 1 .nA. 1 n 
l.Uo. 1U 


M/^r 17 77 

JNLjL.1 000 


None 




1 s7 

loo 




2.2 


"7A 
/O 


07 .77. 17 

Uo.oz. 1 /.o 


7n. ACi~ A A 

oU. 4y.4o 


None 


No SMM/MM source 




1 O 1 A 

lyio 


1 "7 
1 . / 


77 

/ / 


07 .71.71 O 

Uo.ol .21.U 


7n./1 C7Q 

oU.4o.2o 


None 


m77Q7 1 7n7C 

1Uo2o2+oUoj 


Dnnlilllai- at r,1 /1 OO/I \- A/I^.tta JP, A «^»-a /OOO 1 \ 

riacniJJer et at. { iyv4j, Motte oZ Anaie (2UU1 ) 


0O7 

oy2 


1 1 
1 . 1 


78 


03:47:41.6 


32:51:44 


None 3 


B5 IRS1 


Motte & Andre (2001) 


424 


7 


79 


03:47:39.1 


32:52:18 


None 5 


None 


335 


? 


80 


03:26:37.6 


30:15:24 


None 


No SMM/MM source 




302 


0.2 


81 


03:43:50.8 


32:3 :24 


None 


No SMM/MM source 




253 


0.9 


82 


03:33:13.1 


31:19:51 


None 


None 




219 


2.4 


83 


03:32:48.9 


31:09:40 


None 


None 




188 


0.3 



Table A.l: continued. 



1NO. 


D A 

KA 2000 

hh:mm:ss 


^^2000 

dd:mm:ss 


L- luster ' 


Name 


Reference^' 


r S50 

mjy/bm 


i\ Km s 


O A 

54 




J1:U4:dd 


xTau .-.6 

None 


None 




1 OA 

184 


U.z 


Q C 

53 




^ 1 .nn. c q 
J l.UU. J J 


None 


Uj/j4+3UjU 


rv, ,t , a( - „i /100Q\ 

Dent et al. ( iyyo ) 


1 / / 


1 .J 


86 


03:26:30.9 


30:32:28 


None 


None 


176 


1.5 


87 


03:32:21.5 


31:05:08 


None 6 


None 




166 


0.4 


88 


03:31:31.6 


30:43:32 


None 


None 




165 


2.2 


89 


03:32:25.9 


30:59:05 


None 


None 




164 


2.2 


90 


03:45:16.5 


32:04:47 


None 


None 




163 


0.8 


91 


03:29:23.3 


31:36:08 


None 


None 




154 


? 



(1) If the core satisfies the cluster membership criterion (Sect.[3) then the cluster name is given here. 

<2) References are to the earliest previous detection of the submillimetre source. Full cross-identifications with sources at other wavelengths are 
non-trivial due to differing beam sizes and are left for a future analysis, but to indicate an obvious likely counterpart (e.g., IRAS), 'No mm/smm 
source' is marked in the 'Name' column instead of 'None'. 



(3) SCUBA 850 fan peak flux (14" beam) 

<4) £18q integrated intensity (1' beam averaged) at the 850 yum peak. A question mark in this column indicates a SCUBA peak which lies outside 

the C ls O map. 

<5) B5 pair 

<6) B 1 west pair 



